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DTMF Tone Generation and Detection 


ABSTRACT 
This application report describes the implementation of a dual-tone multiple 
frequency (DTMF) tone generator and detector for the TMS320C54x. The author 
provides some theoretical background on the algorithms used for tone generation 
and detection, and documents the actual implementation in detail. Finally, code is 
benchmarked in terms of speed and memory requirements. 


Introduction 


A DTMF codec incorporates an encoder that translates key strokes or digit 
information into dual-tone signals, as well as a decoder that detects the 
presence and the information content of incoming DTMF tone signals. Each 
key on the keypad is identified uniquely by its row frequency and its column 
frequency (see Figure 1). The DTMF generating and decoding scheme is 
not very computationally extensive and can be handled easily by a DSP 
concurrently with other tasks. This report describes an implementation of 
the DTMF codec on the TMS320C54x, a fixed-point DSP designed by Texas 
Instruments (TI™)t especially for telecommunication applications. 


Colin riegueney Group 


Normal Key- 1209 Hz 1336 Hz 1477 Hz 1633 Hz 
pad 
Row- 
Frequency 
Group 


Figure 1. Touch-Tone Telephone Keypad 
(A row tone and column tone are associated with each digit.) 


T Tl is a trademark of Texas Instruments Incorporated. 


DTMF Tone Generation and Detection 


DTMF Tone Generator 


2 DTMF Tone Generator 


2 


The encoder portion and tone generation part of a DTMF codec are based 
on two programmable, second-order digital sinusoidal oscillators, one for 
the row tone and one for the column tone. Two oscillators, instead of eight, 
facilitate the code and reduce the code size. Of course, for each digit that 
is to be encoded, each of the two oscillators needs to be loaded with the 
appropriate coefficient and initial conditions before oscillation can be 
initiated. Since typical DTMF frequencies range from approximately 700 Hz 
to 1700 Hz, asampling rate of 8 kHz for this implementation is within a safe 
area of the Nyquist criteria. Table 1 specifies the coefficients and initial 
conditions necessary to generate the DTMF tones. Figure 2 displays the 
block diagram of the digital oscillator pair. 


For more detail, see Appendix A, which provides some theoretical 
background and a guideline for determining coefficients and_ initial 
conditions for digital sinusoidal oscillators. 


Tone duration specifications by AT&T state the following: 10 digits/sec is the 
maximum data rate for touch-tone signals. For a 100-msec time slot, the 
duration for the actual tone is at least 45 msec and not longer than 55 msec. 
The tone generator must be quiet during the remainder of the 100-msec time 
slot. 


Table 1. Coefficients and Initial Conditions for Sinusoidal Oscillators 


0.82263 


0.58206 
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Figure 2. Two Second-Order Digital Sinusoidal Oscillators 
(Program-flow description of the DTMF generator) 


For the following description of the program flow, it is helpful to consult the 
flowchart shown in Figure 3. Essentially, the series of keypad entries (digits) 
are translated into a series of dual-tones of certain duration that are 
interrupted by pauses of certain duration. Later, the dual-tones enable the 
decoder to identify the associated digits. But, the pauses also are necessary 
to discriminate between two or more identical digits entered successively. 


The program flow, therefore, incorporates two tasks that are swapped after 
certain time intervals. One task (the “tone task”) generates dual-tone 
samples and the other (the “quiet task”) generates pause samples. Each 
task is assigned a certain duration that is controlled by a timer variable. At 
the end of each task, the task has to initialize the timer variable and the 
task-name (tone or quiet) for the next task to be invoked. At the end of the 
quiet task, one very important component is added: A new digit is retrieved 
from the digit buffer and is unpacked. Unpacking means that the digit is 
mapped to the row/column tone properties (oscillator coefficients, initial 
conditions) and pointers are loaded, pointing to the appropriate locations in 
the oscillator property table. 
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The entire program flow is synchronized to the receive-interrupt service 
routine, which provides a perfect clock for real-time processing and constant 
sample output. On completion of the RINT_ISR, the task scheduler is 
invoked, which determines the particular task (tone or quiet) that needs to 
be executed. Both tone task and quiet task check on the timer variable to 
determine if the end of the task duration is already reached. If not, a tone or 
quiet sample, respectively, is generated. If the end of the task duration is 
reached, the next task name and duration is initialized and starts to execute 
with the completion of the next RINT_ISR. The quiet task, additionally, 
unpacks the next digit at the end of its duration. 


Task Scheduler “Quiet Task” “Tone Task” 


Set TASK = 1 
TIMER = 
TIME_LG 


Set TASK =0 
TIMER = 
TIME_SH 


Decrement Decrement 
TIMER TIMER 


Load New Digit 
Increment 
DIGIT_PTR 


Generate 
Dual-Tone 
Sample 


Unpack New 
Digit 


Generate Quiet 
Sample 


Figure 3. Flowcharts of the DTMF Encoder Implementation 
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3 DTMF Tone Detector 


The task to detect DTMF tones in an incoming signal and to convert them 
into actual digits is certainly more complex than the encoding process. The 
decoding process is by its nature a continuous process, meaning it needs 
to continually search an incoming data stream for the presence of DTMF 
tones. 


3.1. Collecting Spectral Information 


The Goertzel algorithm is the basis of the DTMF detector. This method is a 
very effective and fast way to extract spectral information from an input 
signal. This algorithm essentially utilizes two-pole IIR type filters to compute 
DFT values effectively. It is, thereby, a recursive structure (always operating 
on one incoming sample at a time), as compared to the DFT (or FFT) which 
needs a block of data before being able to start processing. The IIR structure 
for the Goertzel filter incorporates two complex-conjugate poles and 
facilitates the computation of the difference equation by having only one real 
coefficient. For the actual tone detection, the magnitude (here, squared 
magnitude) information of the DFT is sufficient. After a certain number of 
samples N (equivalent to a DFT block size), the Goertzel filter output 
converges towards a pseudo DFT value vk(n), which can then be used to 
determine the squared magnitude. See Figure 4 for a short mathematical 
description of the algorithm. More detail is provided in Appendix B. 
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Goertzel Algorithm in short: 


1. Recursively compute for n = 0 ..N 
v, (n) = 2cos (2 k) - V, (n-1)-V,, (n-2) + x(n) 
where V, (-1) = 0 V, (2) = 0 


X(n) = input 


2. Compute once every N 


IX (kK) = Ye (N) Yi" (N) 


vi2(N) + v2(N-1)-2cos (2 k) v2 (N) v2 (N-1) 


Figure 4. Short Mathematical Description of Goertzel Algorithm 


The Goertzel algorithm is much faster than a true FFT, as only few of the set 
of spectral line values are needed and only for those values are filters 
provided. Squared magnitudes are needed for eight row/column 
frequencies and for their eight-second harmonics. The second harmonics 
information later enables discrimination of DTMF tones from speech or 
music. Table 2 contains a list of frequencies and filter coefficients. The 
choice of Nis mainly driven by the frequency resolution needed, which sets 
a lower boundary. Nalso is chosen so that (k/N)fs most accurately coincides 
with the actual DTMF frequencies (see Table 1) assuming ks are integer 
values and fs is a sampling frequency of 8 ksps. 


Table 2. List of Frequencies and Filter Coefficients 


1st HARMONICS 2nd HARMONICS 
(N= 102) fs =8 ksps (N=102) fs = 8 ksps 


COEFFICIENT FREQUENCY COEFFICIENT 
cos(2pi k/N) (k/N)fs/Hz cos(2pi k/N) 
27860 aoe 18 1412 (+ 1.2% ie — 
784 (+ 1.9%) 26745 | 1540. | 1569 (+ 1.9% | tos 


( 
11 | 863 (+ 1.3%) 25529 1704 1725 (+ 1.3% 7010 
re Pi | 941 (0.0%) | 24216 | 1882 | Ete peat a (+0.0%) | 3023, 
1209 ee 19747 | 2418 | 2431 (+0.5%) | -10891 
[3266 | 42 | 2004 (2 0.9% 
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The graphs of Figure 5 (from a Mathcad simulation) show how the Goertzel 
filters resonate with an input signal of matching row/column frequency 
(case 1) and how they do not resonate when the input signal does not match 
their frequency (case 2). Appendix C contains the Mathcad simulation. 


Squared 858 Hz Squared 800 Hz 
Magnitude Magnitude 1483 Hz 
A A 


1483 Hz 


Figure 5. Mathcad Simulation 
(Case 1: 858 Hz, 1483 Hz — digit = 9) 
(Case 2: 800 Hz, 1483 Hz —> digit = undefined) 


3.2 Modifications to the Goertzel Algorithm 


The coefficients of Table 2 reflect the coefficients needed to recursively 
compute the true DFT. The evenly spaced frequency bins of a true DFT 
present an inherent drawback in the DTMF tone-detection process. The 
DFT frequency bins mostly deviate from the true DTMF frequency by an 
amount in the range of up to 2% off center frequency. To be able to meet the 
acceptable bandwidth specifications, a modification of the algorithm departs 
from the true DFT and tunes frequency bins exactly with the DTMF tone 
frequencies. This modification gives up the DFT property of evenly spaced 
frequency bins; and with that, takes two calculated risks: (1) A frequency bin 
is possibly moved inside of the mainlobe of its neighboring frequency bin. 
Therefore, neighboring frequency bins can affect each other. Note that the 
mainlobe of the continuous magnitude spectrum of a rectangular windowed 
sinewave (window is N wide) is exactly the distance of 2 DFT frequency bins. 
(2) This is especially true when column frequencies and 2nd harmonics of 
row frequencies lie close to one another. Note that column frequencies and 
2nd harmonics of row frequencies share the same frequency band. 


DTMF Tone Generation and Detection 7 


DTMF Tone Detector 


To minimize these risks, two rules were applied: (1) Frequency bins are 
spaced apart by more than 0.9/fs. (2) If rule 1 cannot be achieved, ist 
harmonics frequency bin placement has a higher priority than 2nd harmonic 
bin placement. With these rules applied, the coefficient table given in Table 2 
is modified. All necessary fundamental frequency bins essentially can be 
tuned to their DTMF frequencies. A small frequency deviation is needed to 
be accepted only in two cases of second harmonics. The new method of 
using tuned frequency bins facilitates the checking of acceptable 
bandwidths of DTMF tones and helps meet the acceptable bandwidth 
specifications. The modified table of coefficients is given in Table 3. 


Table 3. List of Tuned Frequencies and Modified Filter Coefficients 


1st HARMONICS 2nd HARMONICS 
(N= 102) fs = 8 ksps (N=102) fs = 8 ksps 
DTMF k FREQUENCY COEFFICIENT 2nd FREQUENCY COEFFICIENT 
(k/N)fs/Hz cos(2pi k/N) 


(k/N)fs/Hz cos(2pi k/N) harm. 


697 27980 1394 | 17.93 | 1406 (+ 0.9%) 14739 


770 26956 1540 | 19.72 | 1546 (+ 0.4%) 11414 
852 25701 1704 | 21.72 1704 7549 


2418 —10565 
1336 16325 2672 | 34.07 2672 —16503 


1477 13085 2954 | 37.66 2954 —22318 
1633 | 20.82 1633 9315 3266 | 41.64 3266 —27472 


3.3 Validity Checks 


Once the spectral information (in the form of squared magnitude at each of 
the row and column frequencies and their second harmonics) is collected, 
a series of tests need to be executed to determine the validity of tone and 
digit results. 


A first check makes sure the signal strength of the possible DTMF tone pair 
is sufficient. The sum of the squared magnitudes of the peak spectral row 
component and the peak spectral column component needs to be above a 
certain threshold (THR_SIG). Since already small twists (row and column 
tone strength are not equal) result in significant row and column peak 
differences, the sum of row and column peak provides a better parameter 
for signal strength than separate row and column checks. Tone twists are 
investigated in a separate check to make sure the twist ratio specifications 
are met. 
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The spectral information can reflect the types of twists. The more likely one, 
called “reverse twist”, assumes the row peak to be larger than the column 
peak. Row frequencies (lower frequency band) are typically less attenuated 
as than column frequencies (higher frequency band), assuming a low-pass 
filter type telephone line. The decoder, therefore, computes a reverse twist 
ratio and sets a threshold (THR_TWIREV) of 8 dB acceptable reverse twist. 
The other twist, called “standard twist’, occurs when the row peak is smaller 
than the column peak. Similarly, a “standard twist ratio” is computed and its 
threshold (THR_TWISTD) is set to 4 dB acceptable standard twist. 


The program makes a comparison of spectral components within the row 
group as well as within the column group. The strongest component must 
stand out (in terms of squared amplitude) from its proximity tones within its 
group by more than a certain threshold ratio (THR ROWREL, 
THR_COLREL). 


Finally, the program checks on the strength of the second harmonics in order 
to be able to discriminate DTMF tones from possible speech or music. It is 
assumed that the DTMF generator generates tones only on the fundamental 
frequency; however, speech will always have significant even-order 
harmonics added to its fundamental frequency component. This second 
harmonics check, therefore, makes sure that the ratio of the second 
harmonics component and the fundamental frequency component is below 
a certain threshold (THR_ROW2nd, THRCOL2ndq). If the DTMF signal pair 
passes all these checks, we say a valid DTMF tone pair, which corresponds 
to a digit, is present. 


We now need to determine if the valid DTMF tone information contains 
stable digit information. This is done by mapping the tone-pair to its 
corresponding digit and comparing it with the previously detected digit. We 
call the digit information stable if it has been detected twice successively. 


Finally, we compare the detected digit with the previous-to-last digit. Only if 
the last digit was preceded by a pause do we accept the current digit as a 
valid digit. The detector is then forced into a state where it waits for a pause 
before being able to accept a new digit. This last step is necessary to ensure 
the discrimination of identical keystrokes succeeding one another. 
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3.4 Program Flow Description of the DTMF Detector 
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This DTMF decoder implementation uses a task scheduler in a similar 
fashion as the encoder, again synchronized to the completion of the receive 
interrupt service routine (ISR). The receive ISR is continuously filling a buffer 
of N = 102 words with incoming data. Once the buffer is full, it signals the full 
buffer by setting the input data status word, indatastat, to 1. The task 
scheduler essentially keeps polling indatastat; and once it recognizes a full 
input data buffer (indatastat = 1), it initiates the DTMF detection process by 
calling the appropriate functions. (Consult the flowcharts in Figure 6 and 
Figure 7 might be very helpful at this point.) 


First, the input data status word is reset to signal that the detection process 
has been initiated. Then, the content of the input data buffer is copied into 
an intermediate buffer for processing (double buffering). All the detection 
functions will then operate on the intermediate buffer. The gain control 
function attenuates strong signal inputs and protects the succeeding 
functions from the overflow of the accumulators. Next, the Goertzel filters 
are executed. Since the preceding gain control ensures that overflow cannot 
occur, overflow checking is removed and optimized loops allow fast 
execution. The output of the Goertzel function are the delay states of the 
16 filters which are collected in an array. 


On completion of the Goertzel function, the digit validation checks are 
invoked. The spectral information is computed from the filter delay states 
and collected in an energy template. For the next round of execution, the 
filter delay states are initialized to zero. The energy template is then 
searched for row and column energy peaks. From then on, the detector 
essentially operates in two modes: the tone/digit detection mode or the 
pause detection mode. In the tone/digit detection mode, the detector 
searches for DTMF tone presence and executes all the digit validation tests. 
In the pause detection mode, DTMF tone detection is disabled and the 
decoder first has to await a pause signal. Tone/digit or pause modes are 
controlled by the detectstat variable. Digit validation checks include signal 
strength, reverse and standard twist, relative peaks, second harmonics and 
digit stablility. With the successful completion of these tests, the valid digit 
is stored into the digit output buffer. 
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Buffer Ready 


Indatastat = = 
? 


Yes 
Set 
Indatastat = 0 
Gain Control Attenuate 
Strong Signals 
Goerizel Filters 
for 1st and 2nd Harmonics 


Digital Validation Checks 


Figure 6. Flowcharts of the DTMF Decoder Implementation 
(Part 1 of 2) 
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Task #2 CG) 
Energy Computation 
Initialize filter taps to zero 


Search for energy maxima 


Detector status 
DTMF_STST 
==1? 


Tone/Digit 


D ion M 
Pause Detection Mode election Mods 


Is this a 
Pause 
? 


set 


DTMF_STAT = 


Valid Tones! 


Map detected tone pair to digit 
Store digit in buffer 


Same 
digit as last 
? 


Valid Digit! 


Store digit as valid result 


Figure 7. Flowcharts of the DTMF Decoder Implementation 
(Part 2 of 2) 
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3.5 Multichannel DTMF Tone Detection 


The software is written as C-callable reentrant functions. This enables the 
user to set up multichannel tone detectors in C without adding significant 
additional code. In order to facilitate and structure multichannel applications, 
the code uses a structure to hold all the global variables and pointers to 
various arrays for a single channel. All the user has to do is to define a 
structure for each channel and initialize it properly. A pointer to the structure 
of a desired channel is passed to the C-callable functions and the detection 
takes place. 


DTMF Tone Generation and Detection 13 


Speed and Memory Requirements 


4 Speed and Memory Requirements 


Table 4 and Table 5 summarize the speed and memory requirements of the 
DTMF encoder/decoder implementation. The encoder, as well as the 
decoder, occupies only a very small portion of the computing capability of 
this 50-MIPS DSP. The maximum MIPS count for the actual DTMF encoder 
is approximately 0.45 MIPS. The isolated DTMF decoder uses 
approximately 0.82 MIPS. These speed specifications include all 
processing necessary after the completion of the receive interrupt service 
routine. Also, a sampling rate of 8000 samples/sec is assumed. The MIPS 
performance of the DTMF decoder is achieved using a buffered concept 
instead of a sample-by-sample concept. The speed-critical Goertzel-DFT 
function has the 16 filters coded inline for fast execution. Due to an improved 
gain control, the typical overflow check within the goertzel routine was 
omitted. A drawback of the buffered concept is the need for additional data 
memory due to double buffering of data. The amount of data memory is 
given, to a large degree, by the size of the input data buffers, which are set 
to N = 102 length in this application. Table 4 and Table 5 summarize the 
benchmarks for the DTMF encoder and decoder. 


Table 4. Speed and Memory Requirements for the DTMF Encoder 


MAX 
PROGRAM DATA CYCLES 
MODULE NAME TASKS INCLUDED MEMORY MEMORY PER 
SAMPLE 
DTMF Encoder Task Scheduler 
(sample by sample) pene tae 
pre RY p Quiet Task 
HW/SW Initialization de 
V0 swinit 
RINT_ —e 


TOTAL 


DATA MAX 
MODULE NAME C FUNCTIONS phocnas MEMORY CYCLES 
MEMORY | CHANNELS) | PER 102 

SAMPLES 


some overhead 
DTMF Decoder Gain Control 
(buffered) Goertzel-DFT 

DTMFChecks 


HW/SW Initialization Sig 
inits 


vO C-environment 


TOTAL pO 
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5 Performance 


Two well-known tests to evaluate the performance of DTMF decoders are 
available through MITEL and Bell Communications Research (Bellcore), 
both of which supply the associated test tapes and test procedures. 


5.1 MITEL Tests 


The DTMF tone decoder has been tested using the MITEL test procedures. 
For a preliminary test, a set of files with the digitized signal data for the 
various MITEL tests was acquired through TI internal sources. The files 
essentially contain a subset of the signal data of the real MITEL test tapes. 
Since none of the files contained more than 32K words of data, the files were 
reformatted as “.inc” files and then included in the source code. At link time, 
the data of a given test file was mapped into the external memory of the 
TMS320C54xEVM, occupying at maximum 32K words of data space. The 
code was then tested in real-time using the contents of the 32K-word test 
buffer as its input source. During this preliminary test, the various decoder 
thresholds were set to proper levels. 


The complete MITEL test was then executed according to the test 
procedures specified in the MITEL test document. A digital audio tape (DAT) 
recording of the test tapes was used as input signal source. The decoder 
was executed on the TMS320C54xEVM utilizing a TMS3820AC01 analog 
interface to convert the incoming signal into the digital domain. The MITEL 
test essentially has two sections. The first section measures the DTMF tone 
decoder in terms of recognition bandwidths (RBW), recognition center 
frequency offset (RCFO), standard twist, reverse twist, dynamic range (DR), 
guard time and signal-to-noise ratio (SNR). The second section, called 
“talk-off test”, consists of recordings of conversations on telephone trunks 
made over a long period of time and condensed into a 30-minute period. In 
this section, the decoder’s capabilities to reject other sources such as 
speech and music is measured. MITEL specifies a maximum of 
30 responses of the DTMF decoder as acceptable speech-rejection. 
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Table 6 summarizes the MITEL test results. With the exception of the 
recognition bandwidth results for the low-band frequencies, the DTMF 
decoder passes all the tests and exceeds the specifications. The slightly 
higher relative bandwidths for the low-band frequencies were accepted in 
favor of a higher robustness. Essentially, the signal strength threshold 
THR_SIG2 determines to a large part the recognition bandwidths. 
THR_SIG2 has been setto leave slightly more room for digit acceptance and 
higher robustness. The threshold settings for acceptable twists helped pass 
the twist test and exceed the specifications. The dynamic range of the 
decoder of 28 dB is better than the specification. The guard time of 30 ms 
leaves enough head room to safely detect digits with a specified 45-ms tone 
duration. The decoder was able to correctly detect all 1000 tone bursts for 
each of the given noise environments of —24 dBV, —18 dBV and —12 dBV 
AWGN. When the decoder was exposed to the 30 minutes of speech and 
music samples, it did not respond a single time and, thereby, exceeded the 
MITEL talk-off specification of 30 permissible responses by a significant 
margin. 
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Table 6. MITEL Test Results 


BW TESTS FREQUENCY RBW% 
specification Po 1.5% < RBW < 3.5% 


697 Hz 5.8% 


Low band 
High band 
TWIST TESTS REV TWIST 
specification >8dB 
DIGIT 16 8.5 dB 
DATESTS 
specification 
DIGIT 
DIaIrs 
Diair9 
oIair 16 
GUARD TIME MIN TONE TIME 
specification Pt 45 ms 


SNR TESTS RESULT 


| NOISE 
specification 
specification 
result very robust 
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5.2 BELLCORE Talk-Off Test 


Through TI internal sources, the Bellcore series-1 Digit Simulation Test 
Tapes for DTMF receivers were available for testing. These tapes consist 
of six half-hour sequences of speech samples, designated part 1 through 
part 6, which are known to contain energy at or near valid DTMF frequency 
pairs. This test exhaustively measures the speech-rejection capabilities of 
DTMF receivers in telecommunication systems. There are over 
50,000 speech samples, some of which are music samples. It is estimated 
that the six parts of the series-1 tapes are equivalent to the exposure of one 
million customer dialing attempts in a local central office. In other words, 
exposing a DTMF receiver to all the speech samples in series-1, will produce 
the same number of digit simulations that the receiver would experience if 
it were exposed to customer speech and room noise present during network 
control signaling on one million calls. The Bellcore talk-off test is far more 
exhaustive than the MITEL talk-off test. 


The test setup was identical to the one used for the MITEL talk-off. Table 7 
summarizes the test results for the Bellcore talk-off. In the three hours of 
testing, the decoder responded to digit simulations only in six cases. This is 
far less than the specifications required to pass the talk-off. The decoder 
proved to be very robust in terms of its soeech-rejection capabilities. 


Table 7. Bellcore Talk-off Test Results 
SPECIFICATION RESULTS 


per 1,000,000 calls 

i ' oes : . ite pe = 333 FESPONses . aia 
‘ per 1,000,000 calls 

isi ' alii i . ne Ss 500 FESPOnses ° aa 
per 1,000,000 calls 

parts 1 through 6 (3 hours) 0-9, *, #, A, B, C, D < 666 responses 6 responses 
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6 Summary 


DTMF tone encoding and decoding concepts and algorithms were 
described in detail. Further theoretical background is provided in the 
appendix. The DTMF encoder and decoder implementations were 
explained, and the associated speed and memory requirements were 
presented. The DTMF tone decoder has been tested according to the 
MITEL and BELLCORE test specifications and the results are documented. 
It is important to note that the decoder was implemented as reentrant, 
C-callable functions, which facilitate setting up a multi-channel DTMF 
decoder system. The code is modular and easy to integrate into any given 
telephony application. The decoder algorithm was optimized to meet the test 
specifications as well as offer a very attractive MIPS count far below 1 MIPS 
per channel. 
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Background: Digital Oscillators 


Appendix A_ Background: Digital Oscillators 


A.1 Digital Sinusoidal Oscillators 


A digital sinusoidal oscillator can, in general, be viewed as a form of a 
two-pole resonator for which the complex-conjugate poles lie on the unit 
circle. It can be shown that the poles of a second-order system with system 
function 


bo 


a (A.1) 
1+ a2! + aoz* 


H(z) 


with parameters 


by = Asinwy 
2 COS Wy (A.2) 


a) 
Il 


a, = 1 
are exactly located at the unit circle. That is 

Py 5 = e*/% (A.3) 
The discrete-time impulse response 

h(n) = Asin ((n + 1) wo) - u(n) (A.4) 


corresponding to the above second-order system clearly indicates a clean 
sinusoidal output due to a given impulse input. Therefore, this system can 
be termed a digital sinusoidal oscillator or digital sinusoidal generator. 


For the actual implementation of a digital sinusoidal oscillator, the 
corresponding difference equation is the essential system descriptor, given 
by 


y(n) = -ayy(n-1)-asy(n-2) + b)d(n) (A.5) 


where initial conditions y(—1) and y(—2) are zero. Note that the impulse 
applied at the system input serves the purpose of beginning the sinusoidal 
oscillation. Thereafter, the oscillation is self-sustaining, as the system has 
no damping and is exactly marginally stable. Instead of applying a delta 
impulse at the input, let the initial condition y(—2) be the systems oscillation 
initiator and remove the input. With this in mind, the final difference equation 
is given by 


y(n) = 2c0s ay: y(n—-1)-y(n-2) (A.6) 
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where 
Val}; =0 
y(-2) = -Asinay (A.7) 
Wg = 2nfo/fs 


with f being the sampling frequency, fo being the frequency and A being the 
amplitude of the sinusoid to be generated. Note that the initial condition y(—2) 
solely determines the actual amplitude of the sinewave. 
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Appendix B_ Background: Goertzel Algorithm 


Goertzel Algorithm 


As the first stage in the tone-detection process, the Goertzel algorithm is one 
of the standard schemes used to extract the necessary spectral information 
from an input signal. Essentially, the Goertzel algorithm is a very fast way 
to compute DFT values under certain conditions. It takes advantage of two 
facts: 


1. The periodicity of phase factors {WK} allows the expression of the 
computation of the DFT as a linear filter operation utilizing recursive 
difference equations. 


2. Only a few of the spectral values of an actual DFT are needed (in this 
application, there are eight row/column tones plus an additional eight 
tones or corresponding 2nd harmonics). 


Keeping in mind that a DFT of size Nis defined as 


N-A1 
X(k) = 5° x(m)e Fem (B.1) 


m=0 


it is possible to find the sequence of a one-pole resonator 
N-1 ™ 
y,(n) = >» x(m) e/'n k(n) (B.2) 
m=0 


which has a sample value at n= N coinciding exactly with the actual DFT 
value. In other words, each DFT value X(k) can be expressed in terms of the 
sample value at n= Nresulting from a linear filter process (one-pole filter). 


It can be verified that 


= x(m) elWhN 9 apkm (B.3) 


ll 
= 
3 
© 
at 

2| 
a> 
3 


The difference equation corresponding to the above one-pole resonator 
sequence (B.2), which is essential for the actual implementation, is given by 


DTMF Tone Generation and Detection 


B-1 


Background: Goertzel Algorithm 


B-2 


(B.4) 


with y(-1) = 0 and pole location p = e7w* . Being a one-pole filter, this 
recursive filter description yet contains complex multiplications, not very 
convenient for a DSP implementation. Instead, by using atwo-pole filter with 


; 2m 
; 2n es ae 
complex conjugate poles p; 5 = e +/'v" and only real multiplications in its 
difference equation, 


N 
where y%x(-1) and y%.(—2) are zero. 


v, (m) = 2cos (22 k) vy, (n-1)-V, (7-2) + x(n) (B.5) 


In the Nth iteration, only a complex multiplication is needed to compute the 
DFT value, which is 


X(k) = y,(N) = v,(N)-e Fk v, (NA) BS) 


However, the DTMF tone-detection process does not need the phase 
information of the DFT; squared magnitudes of the computed DFT values, 
in general, suffices. After some mathematical manipulation, it is found that 


= VAN) Y,," (N) 


= v,2(N) + v,2 (N-1)-2cos (2 k) v,.2(N) v,2(N-1) 


|X(k)| 
(B.7) 


In short: For the actual DSP implementation, equations (B.5) and (B.7) are 
used to retrieve the spectral information from the input signal x(n) for further 
evaluation. Note that equation (B.5) is the actual recursive linear filter 
expression, which is looped through for n = 0 .. N. Equation (B.7) is 
computed only once every N samples to determine the squared 
magnitudes. 
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Appendix C DTMF Tone Generator Code Listing 


DTMF Tone Generator Executable on a TMS320C54x EVM 


Kak kkk KR ARARK KKK KKK KKK KKK kk kkk kkk kkk kkk kkk kkk kkk kkk Kk KKK KKK KKK KKK 


INC. 1996 * 


KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KK 


08/20/96 


KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKKKKKK KKK KK KKK KKK 


DTMF tone generator - 
dtmf100e.asm 

This file contains the code for a 
DTMF tone generator 

(TMS320C54x EVM version) 


* 
* 
* 
* 
* 
* 
* 
* 
* 
* 


; define global vars and consts 
; Gefine global tables and buffers 


* (C) COPYRIGHT TEXAS INSTRUMENTS, 

* Program Name: 

2 File Name: 

* File Description: 

* 

* 

* 

* Author: Gunter Schmer 

* Date: 08/20/96 

- Revision: 2.0 

zx Latest working date: 

* 
-title “"dtmf encoder” 
-mmregs 
-include “"varse.inc” 
-include “"sectse.inc” 
-include "globalse.inc” 


; globalized labels 


Kkk kk RARKRK KKK KKK KKK KKK KK KKK KK KKK KK KK KK Kk KK 


* CODE STARTS HERE 


KKK KKK KKK KKK KK KKK KKK KKK KKKKKK KKK KK KKK KKK KK 


-text 
START 

ssbx 

ld 


call 


ssbx 
rsbx 
ssbx 


stm 
stm 
rsbx 
ld 
mvmd 
mvdm 


INTM 
#0,DP 


hwinit 


SXM 
OVM 
FRCT 


#0040h, IMR 
#00c8h, IFR 
INTM 
#aic_conf,DP 
DRR1, rcv 
tra,DXR1 


’ 


global interrupt disable 

initialize data pointer 

Initialize Hardware 

hardware initialization subroutine 
Initialize Software 

data sign ext. before usage 

no saturation of accu if overflow 
fractional mode bit, left shift of 
multiplier to compensate for extra sign bit 
enable RINT1 interrupt 

clear all pending serial port or timer interrupts 
global interrupt enable 

DP to variables 

dummy read 

dummy write 


KKK KKK KKK KKK KKK KKK KKK KKK KKK KK KKK KKK KKK KK KK 


a3 Main Program 
Kkkh kkk ARARKRK KKK KKK KKK KKK KKK KKK KK KK KKK KKK KKK 


Main 

st 
st 
st 
stm 
idle 
call 
nop 


next 


#PH_NBR, DIGIT_PTR 


#0, TASK 
#0, TIMER 
#DATA, AR5 
1 

tasks 


next 


; initialize pointer to PHNBR table 
; initialize beginning task 

; initialize timer to zero 

; ARS points to testbuffer 


; actual processing following RINT interrupt 


7 store to test buffer 


; stm #7f£ffh,BK 
p> lid tra,A 
Pst A, *AR5+% 
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done b done 

KKKKK KKK KKK KKK KKK KKK KKK KKK KK KKKKKKKKKKK KKK KKK KKK KKK KKK KKK KKK 
* Task Scheduler: 

* Run this task scheduler section at the beginning of 


* each RINT_ISR process 
* 


* 


KKKKKK KKK KKK KKK KKK KKK KKK KKKKKKKKKKKKKKK KKK KKK KKK KKK KKK KKK KKK 


tasks  cmpm TASK, #00h 
be task1,NTC ; if (TASK!=0) branch to taskl 
; else branch to task0 
taskO cmpm TIMER, #00h 


be task01,NTC ; branch if timer not zero 
mvdm DIGIT_PTR,AR1 ; AR1 points to digit in PHNBR table 
st #01h, TASK ; TASK of next RINT_ISR is taskl 
st #TIME_LG,TIMER ; set timer to long duration 
ld *ARI+,A ; load A with digit 
be done,alt ; done if digit is -1 ! 
mvmd AR1,DIGIT_PTR ; save new pointer to digit in PHNBR table 
call unpack ; unpack digit: A --> T1_OFS, T2_OFS 
b task02 
task0l ld TIMER, A 
sub #1,A 
stl A, TIMER ; decrement timer 
task0O2 call quiet 
b task3 ; branch to task3 
taskl cmpm TIMER, #00h 
be task11,NTC ; branch if timer not zero 
st #00h, TASK ; TASK of next RINT_ISR is task0O 
st #TIME_SH, TIMER ; set timer to short duration 
b task12 
task11 ld TIMER,A 
sub #1,A 
stl A, TIMER ; decrement timer 
task12 mvdm T1_OFS,AR2 ; AR2 is offset for row-tone 
mvdm T2_OFS,AR3 ; AR3 is offset for column-tone 
cali tone 


task3 ret 


KKKKKKKKKKKKKKKKKKKKKKKKKKKKK KKK KKK KK KK KKK 


* Interrupt Service Routines 
KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK KKK KKK KKK KK 
RINT1_ISR: 

ld #rcv, DP 

mvmd DRR1, rcv 

andm #O0fffch,tra 

mvdm tra,DXR1 

popm STi 

popm STO 

rete 

.end 


;varse.inc 
KK KK KKK KKK KKK KKK KKK KKK KKK KK KKK KKK KKKKEKK 


aee Constants 
KKEKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK KK 

TIME_SH .set 440 7 pause duration (55 msec) 
TIME_LG .set 360 ; tone duration (45 msec) 


DAC_OFS .set 000h ; DAC offset 


KKKKKKKKKK KKK KKKKKK KKK KKKKKKKKKKKKKKKKK 


ee Variables 
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Kkk kkk ARARKARK KKK KKK KKK KKK KKK KKK KK KKK Kk KKK 


.bss aic_conf,1 

-bss rev,1l ; veceive variable 

.bss tra,l ; transmit variable 

-bss DIGIT,1 ; current digit 

-bss DIGIT_PTR,1 ; points to current digit in PHNBR table 

-bss TIMER, 1 ; timer counter 

-bss TL ORS. ; offset for tone 1 (row tone) 

-bss T2_OFS,1 ; offset for tone 2 (column tone) 

-bss TASK, 1 ; holds the nbr of the task that is to be executed 


;sectse.inc 
Kaka KKa KKK KKK KK KKK KKK KKK KKK KK KKK KKK KKK KKK 


*x*x Tables 
Kaa K aK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK 


Kkk kk KK KR ARKRK KKK KKK KKK KKK Kk kkk kkk kkk kkk kkk kk kkk kkk kk kkk KK 


* The following table assembles the coefficients and 


* initial conditions for the difference equations of * 
* the digital sinusoidal oscillators. es 
* * 
* In general: - 
* DEQ: y(n) = 2*cos(2pi*f/fs)*y(n-1) - y(n-2) * 
% EC. y(-1) = 0 oe 
* y(-2) = -A*sin(2pi*f/fs) * 
* where A = desired amplitude of sine wave * 
x f = desired frequency of sine wave * 
bas fs = sampling frequency - 
* * 
* Example: x 
* -word cos (2pi*f/fs)*32768 ;coefficient x 
a -word 0 7y(-1) * 
“ -word -A*sin(2pi*f/fs) *32768 7y (-2) - 
* * 
KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KEKE KEKE KKK KKK KKK KKK KKK KKK KKK 

-sect "tbhl_tone” 
TONES -word 27980 7; row coef 

-word 0 7 y(n-1) 

-word -1024*5204/10000 7; y(n-2) 

-word 26956 7; vow 2 coef 

-word 0 7 y(n-1) 

-word -1024*5686/10000 7; y(n-2) 

-word 25701 7 row 3 coef 

-word 0 7 y(n-1) 

-word -1024*6203/10000 7; y(n-2) 

-word 24219 7 row 4 coef 

-word 0 7 y(n-1) 

-word -1024*6736/10000 7; y(n-2) 

-word 19073 ; col 1lcoef 

-word 0 7 y(n-1) 

-word -1024*8132/10000 7; y(n-2) 

-word 16325 7 col 2 coef 

-word 0 7 y(n-1) 

-word -1024*8671/10000 7; y(n-2) 

-word 13085 ; col 3 coef 

-word 0 7 y(n-1) 

-word -1024*9168/10000 7; y(n-2) 

-word 9315 ; col 4 coef 

-word 0 7 y(n-1) 

-word -1024*9587/10000 7; y(n-2) 
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KKKK KKK KKK KKK KKK KKK KK KKK KKK KKK KKK KKKKKKKKKKKKKK KKK KKK KKK KKK KKK KKK 


The following table contains offsets into the 


+ + F F 


KEYS 


-sect 
-word 
-word 
-word 
-word 
-word 
-word 
-word 
-word 
-word 
-word 
-word 
-word 
-word 
-word 
-word 
-word 


tone table tbl_tone. 
The offset for the row-tone is the 
The offset for the column-tone is the 


"tbl_keys” 


0903h ; 
0000h ; 
0003h ; 
0006h . 
0300h ; 
0303h : 
0306h ; 
0600h F 
0603h ; 
0606h H 
0009h ; 
0309h : 
0609h ; 
0909h F 
0900h ; 
0906h ; 


ror 
myer 
ror 
13° 
rar 
mor 
6! 
ares 
rgr 
ror 
ar 
BI 
rc 
pr 
Er 
Pr 


example: 


= leer 


= Hr 


(upper byte) 
(lower byte + #12) 


KKKKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKKKKKKKKKKK KKK KKK KKK KKK KKK KK KK 


0903h --> 3.row, 


of each word 


* 
* 
* 
* 
* 


1.column 


KKEKKKKKKKKK KKK KK KKK KKK KKK KKK KKK KKK KKK KKK KKK KK KKKKEKKKKKKKK KKK KKK HK 


The following table contains the phone number to be 


PH_NBR 


is the foll nbr: 


KKKKKKKK KKK KKK KKK KKK KKK KKK KKK KKKKKKKKKKKKKKK KKK KKK KKK KKK KK KK KKK 


-sect 
-word 
-word 
-word 
-word 
-word 
-word 
-word 
-word 
-word 
-word 
-word 
-word 
-word 
-word 
-word 
-word 
-word 


-word 


427* #A 3 


"tbhl_phnr” 

1,2,3,4,5, 6, 7,8, 9, 
1,2,3)4,57 6,7, 879, 
1,2,3,4,5,6, 7,8, 9, 
1,2,3,47 57 6, 7,87 97 70. 
1,2,374757, 677,879, 70 
1,2,3,4,5,6, 7,8, 9, 70. 
17273747576, 17.849, 70 
1,2,3,475976),7,8)9 
1,2,3, 475,76, 7,8, 9, 70 
1,2,3,4,5,6,7,8,9, 70 
1,2,3,4757 677,879, 
1,2,3,4,5,6, 7,8, 9, 70. 
1,2,3747576, 778,79, 70 
1,2,3,4,5,6,7,8, 9, 
1,2,3,4,5,6, 7,879, 70 
1,2,3,4,5,6,7,8,9 0 


~ 
~ 
~ 
~ 
~ 
~ 
~ 


, 


=1 


encoded. The format is as follows 


tf 


4,2,7,0Eh, OFh, 0OAh, 3,-1 


+ + HF 


* 


(-1 terminates the encoding) * 


KKEKKKKKKKKKK KKK KKK KK KKKKKKKEKKEKKKKK KKK KK KKK 


* TEST BUFFER 


KKEKKKKKKKKKK KKK KKKK KKK KKKKKKKKKKKKKK KKK KKK 


DATA 


-usect 


"testbuf”,8000h 


KKKKKKKKKKKKKKKKKKKKKKKKKKKKKK KKK KKK KKK KKK 


* 


BOS 
TOS 


C-4 


Stack setup 


KKEKKKKKKKKK KK KKK KKKKKKKKKKKKKKKKKK KK KKK KK 


-usect 
-usect 


"stack”,20h : 
"stack”,1 2 
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setup stack 


0, OAh, 0Bh, 0Ch, 0Dh, 0Eh, OFh 
0, OAh, 0Bh, 0Ch, 0Dh, 0Eh, OFh 
0, OAh, 0Bh, 0Ch, 0Dh, 0Eh, OFh 
0, OAh, 0Bh, 0Ch, 0Dh, 0Eh, OFh 
0, 0OAh, 0Bh, 0Ch, 0Dh, OEh, OFh 
0, OAh, 0Bh, 0Ch, 0Dh, 0Eh, OFh 
0, OAh, 0Bh, 0Ch, 0Dh, 0Eh, OFh 
, 0, 0Ah, 0Bh, 0Ch, 0Dh, 0Eh, OFh 
0, OAh, 0Bh, 0Ch, 0Dh, OEh, OFh 
0, OAh, 0Bh, 0Ch, 0Dh, OEh, OFh 
0, OAh, 0Bh, 0Ch, 0Dh, 0Eh, OFh 
0, OAh, 0Bh, 0Ch, 0Dh, OEh, OFh 
0, OAh, 0Bh, 0Ch, 0Dh, OEh, OFh 
0, OAh, 0Bh, 0Ch, 0Dh, 0Eh, OFh 
0, OAh, 0Bh, 0Ch, 0Dh, 0Eh, OFh 
0, OAh, 0Bh, 0Ch, 0Dh, 0Eh, OFh 


Top of stack at reset 
KKEKKKKKKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKKKKKKK KKK KKK 


(*3) 
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* 


(C) 


COPYRIGHT TEXAS INSTRUMENTS, 


INC. 1996 


* 


Kkkk kkk KR ARKRK KKK KKK KKK KKK kkk kkk kkk kkk kkk kkk kkk Kk KKK KKK KKK KK KKK Kk 


Program Name: 
File Name: 
File Description: 


DIMF tone generator 
dtmfsube.asm 
This file contains subroutines for a 


* 


Author: 
Date: 
Revision: 


Latest working date: 
KKK KKK KKK KKK KKK KKK KK KKK KK KKK KKK KKK KKK KKK KKK KKK KKKKKKKKK KK KKK 


DTMF tone generator 
(TMS320C54x EVM version) 


Gunter Schmer 
08/20/96 

2.0 

08/20/96 


-mmregs 


- include 


"globalse.inc” 


* 
* 
* 
* 
* 
* 
* 
* 
* 
* 


KKK KKK KKK KKK KKK KK KKK KKK KKK KKK KKK KKK KKK KKK KKK KEKE KKK KKKKK KKK KKK KK KKK 


* SUBROUTINE: 
* Description: 


* 
* 
% two tones and writes results into variables = 
* T1_OFS, T2_OFS * 
* Uses: AR2 * 
* Input: A 
* Output: none ‘a 
Kkk kk kkk KR ARKRKR KKK KKK KKK KKK kkk kkk kkk kkk kkk kkk kkk kkk kkk kkk kkk kk kk kkk 
unpack stim A, AR2 ; AR2 keys-map offset for unpacking 
ld #DIGIT,DP ; DP to variables 
nop ; pipeline conflict(1) + latency(1) 
ld *AR2 (KEYS),A ; load A with keys-map value 
ld *AR2 (KEYS) ,B ; load B with keys-map value 
and #0£00h,A ; mask out the row portion 
and #000fh,B ; mask out the column portion 
sfitl A,-8 ; vight shift A by 8 
add #12,B ; add 4*3 words offset to point into column portion 
stl A, T1_OFS ; store tone 1 offset (row tone) 
stl B, T2_OFS ; store tone 2 offset (column tone) 
ret 
KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KEKE KKK KEKE KEK KKK KKK KKK KKK KK KKK 
SUBROUTINE: tone * 
* Description: Generates the dual tone samples for DTMF * 
* using offsets T1_OFS, T2_OFS * 
* Uses: AR2, AR3 % 
* Input: none * 
* Output: none * 
KKK KKK KKK KK KKK KKK KKK KKK KKK KKK KKK KKK KKK KEKE KEKE KEK KKK KKKKK KKK KKK KK KK KKK 
tone ld #0,B ; clear B  ----ROW TONE---- 
sub *AR2 (TONES+2),15,B 7 (B) = -(1/2)y(n-2), 
ltd *AR2 (TONES+1) ; load (T) with y(n-1) and y(n-1) 
mpy *AR2 (TONES),A 7 (A) = coef*y(n-1) 
add A,B ; (B) = coef*y(n-1) - (1/2) y(n-2) 
sth B, 1, *AR2 (TONES+1) ; 2*(B) --> y(n-1) 
ld #0,B ; clear B  ----COLUMN TONE---- 
- b tonel ; test only one oscillator 
sub *AR3 (TONES+2),15,B ; (B) = -(1/2)y(n-2), 
ltd *AR3 (TONES+1) ; load (T) with y(n-1) and y(n-1) 
mpy *AR3 (TONES),A 7 (A) = coef*y(n-1) 
add A,B 7; (B) = coef*y(n-1) - (1/2)y(n-2) 
sth B, 1, *AR3 (TONES+1) ; 2*(B) --> y(n-1) 
tonel add *AR2 (TONES+1),15,B ; add two tone samples 
add #DAC_OFS,B ; add DAC offset 
sth B,tra ; write to transmit variable 


unpack 


Maps the key value (in A) 


into two offsets for 


DTMF Tone Generation and Detection 


with AR2 


in high accumulator 
--> y(n-2) 


in high accumulator 
--> y(n-2) 
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ret 
KK KKK KKK KKK KKK KKK KKK KKK KK KKK KKK KKK KKK KKKKKKKKKKKKKKK KKK KKK KK KK KKK 
* SUBROUTINE: quiet x 
* Description: Generates a pause of specified duration * 
* * 
* Uses: AR2, AR3 is 
* Input: TIME_SH i 
* Output: none - 
KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KK KKK KKK KKK KKKKKKKKKKKKKKKKKK KKK KKK 
quiet ld #0,B ; clear B 

add #DAC_OFS,B ; add DAC offset 

sth B,tra ; write to transmit variable 

net 
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Appendix D DTMF Tone Detector Code Listing 
D.1 DTMF Tone Detector Executable on a TMS320C54x EVM 


[KK KK HK KKK I KR IK KA IRA RRA IA RIA IKK RI IK He 


* (C) COPYRIGHT TEXAS INSTRUMENTS, INC. 1996 * 


KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KEK KEKE KEKE KEKE KK KKK KKK KK KKK KKK KK 


* Program Name: DTMF tone decoder 
2 File Name: globals.h x 
x File Description: defines channel information * 
* * 
ia Author: Gunter Schmer % 
= Date: 10/15/96 * 
* Revision: 2.0 * 
* Latest working date: 10/15/96 * 


FR AR RRA A RIAA IRR RA IRR RIA IRA I I / 
/* globals */ 

/* struct used to pass on channel info to functions * 
typedef struct 

{ 


int *indata; /* ptr to indata */ 

int *taps; /* ptr to filter states if 

int *energy; /* ptr to energy template * 

int *digitptr; /* ptr to digit output array * f 

int digitlast; /* last detected digit */ 

int detectstat; /* status of detector */ 

/* 0 - detector waiting for pause af 

/* 1 - detector ready to detect digit 7. 

int err_flags; /* error flags for dtmf checks */ 

/* bO - error signal strength */ 

/* bl - error reverse twist */ 

/* b2 - error standard twist */ 

/* b3 - error row’s relative peak */ 

/* b4 - error col’s relative peak */ 

/* b5 - error row’s 2nd harmonic * f 

/* b6 - error col’s 2nd harmonic yf 

/* b7 - error OVERFLOW of ACCA */ 

/* b8..b15 are zero x] 

} DTMFCHANNEL; 

[RR soe cc ese channel#1 information ------------------------------ kf 
int indatal[102]; /* receive buffer */ 
int indata2[102]; /* process buffer */ 
int digits [128]; /* output buffer for digit results ff 
int taps[32]; /* filter states */ 
int energy[16]; /* energy template */ 
/* used for interrupt servicing and signaling a full buffer xf 
int *indataptr = é&indatal[0]; /* ptr to receive buffer */ 
int indatastat = 0; /* status of receive buffer */ 
/* 0 = buffer not ready to process */ 
/* 1 = buffer ready to process */ 


/* initialized DIMFCHANNEL */ 
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DTMFCHANNEL channell = { é&indata2[0], 
é&taps[0], 
éenergy[0], 
édigits[0], 


OxFFFF, 
1, 
0 }; 

DTMF CHANNEL *ptrchannell = &channell; 
[PP SeesHaaaesaa-= end of channel#1 —---~--~--~---~~---~-------- a A 
/* used for testing */ 
volatile unsigned int *testdataptr = (volatile unsigned int *) 0x8000; 
[ROR KKK KKK KKK KK KKK IKK IK KKK IK KK II IOI IR IK KK 
* (C) COPYRIGHT TEXAS INSTRUMENTS, INC. 1996 * 
KKEKKKKKK KKK KKK KKK KKK KKK KKK KKK KKK KK KKK KKK KKK KKKKKKKKKKKKK KKK KKK 
s Program Name: DTMF tone decoder * 
* File Name: main.c * 
* File Description: calls init functions - 
* calls decoder functions e 
* * 
bs Author: Gunter Schmer " 
* Date: 10/15/96 * 
Revision: 2.0 * 
Latest working date: 10/15/96 * 
ee ee ee eo 
#include “init.h” 
#include “globals.h” 
#include "mmregs.h” 
#include "dtmfsub.h” 
/* prototypes */ 
void tasks (void) ; 
void initArrays (void); 
main () 
{ 

[KOR KK KK AK a OK OK init HW and SW TK AK KK A I  / 

initWSGen (); /* initialize WS generator*/ 

initArrays(); /* initialize arrays * 

initSP1(); /* initialize SP1 */ 

initACO1(); /* initialize AIC * if 

initInts (0x0040) ; /* allow only RINT1 */ 

[RRR KKK KKK KKK KKK synch tasks () ee: RINTO KOR KK KK / 

for (;;) { 
asm(” IDLE Lhe /* wait for interrupts */ 
tasks (); /* run tasks on completion of ISR */ 


} 
void tasks (void) 
{ 
int n; 
if (indatastat==1) { 
/* reset indata status bit */ 
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} 


indatastat = 0; 

/* copy indatal to indata2 */ 
for (n=0;n<102;n++) 

indata2[n] = indatal[n]; 

/* perform gain control */ 


gaincntrl(ptrchannell) ; 


/* perform goertzel */ 


goertzel (ptrchannell); 


/* perform dtmf digit validation 


dtmfchecks (ptrchannell) ; 


void initArrays (void) 


{ 


int n; 

for (n=0;n<102;n++) 
indatal[n] = 0; 
indata2[n] = 0; 


} 

for (n=0;n<128;n++) 
digits[n] = 0; 

for (n=0;n<32;n++) 
taps[n] = 0; 

for (n=0;n<16;n++) 
energy[n] = 0; 


{ 


a 
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KKK KKK KKK KKK KK KKK KKK KKK KKK KKK KEK KKK KKK KEKE KEK KKK KKK KK KKK KKK KKK KK KK 


* 


(C) COPYRIGHT TEXAS INSTRUMENTS, 


INC. 


1996 


* 


Kkhk kk kkk KR KR AR KRK KKK KKK KKK KKK KKK Kk kkk kkk kkk kk KKK Kk Kk KKK KKK KK KKK KKK Kk 


* 


Program Name: 
File Name: 
File Description: 


Author: 
Date: 


Revision: 2.0 


DIMF tone decoder 
dtmfsub.asm 


* 


This file contains the subroutines * 


used for a DIMF tone detector 
(TMS320C54x EVM version) 
Gunter Schmer 
10/15/96 


Latest working date: 10/15/96 


Kkhk kkk kK ARKRKK KKK KKK KK KKK Kh kkk kkk kkk kkk kkk kkk kkk kK KKK KKK KK KKK KKK 


-mmregs 
- include 
- include 


"globals.inc” 
"tables.inc” 


7global labels 
;various tables 


Kkk kk kkk KR AR KRK KKK KKK KKK KK kkk kkk kh kkk kkk kkk kkk kkk kkk kkk kkk Kk KKK KKK KK KKK 


* 


* 


* 


* 


SUBROUTINE: gaincntrl (approx. 705 cycles) 


Input: none 
Output: none 
Uses: 


Description: Scales a block of 


input data in order to ensure 
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no overflow when running the goertzel routine 


approx. 490 cycles per N=102 sample frame 


approx. 5 cycles per sample (negligable) 


gain_power = (1/64) * SUM{0,N-1}[x(n) *x(n)] 


(1/64) * N * signalpower 


gain_lim = (1/64) * (1/N) 

= minimum signalpower (-..dB) 
gain_const = 0.5*sqrt (2) *sqrt (64/N) 
gain_amp = gain_const * sqrt (gain_power) 
gain_lev = 1/N 
gain_scale = gain_lev/gain_amp 
x(n) = x(n) * gain_scale 


* FRCT=1, SXM=1 


KKKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KK KKK KKKKKKK KKK KKK KKK KKK KKK KKK 


-text 
_gaincntrl: 

nop 

pshm STO 

pshm STi 

pshm AR1 

pshm AR6 

pshm AR7 

nop 

frame #-10 

stim A, ARS 

nop 

ld *AR5,A 


stl A, *SP (9) 


ld #0005h,16,A 
or #0505h,A 


sth A, *SP (2) 
stl A, *SP (3) 


ld #00E3h,A 
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;---- local variables ----- 
; *SP(0) = gain_pow (HI) 

; *SP(1) = gain_pow (LO) 

; *SP(2) = gain_lim (HI) 

; *SP(3) = gain_lim (LO) 

; *SP(4) = gain_amp 

; *SP(5) = gain_lev 

; *SP(6) = gain_scale 

; *SP(7) = gain_const 

; *SP(8) = gain_delta 

; *SP(9) = ptr to indata 
jo-- > arguments -=---—--=== 
, (A) = ptr to dtmfchannel 


;value initializations 


7 (A) = 00050000h 

7(A) = 00050505h = 1/(64*N) (Q31) 
7 (A) --> gain_lim (32bit) 

; (AL) = OOE3h = (0.5*sqrt(2))/N 


stl A, *SP (5) ; (AL) --> gain_lev (16bit) 
ld #6564h,A ; (AL) = 6564h = sqrt (64/102) 
stl A, *SP (7) ; (AL) --> gain_const 

ld #4000h,A ; (AL) = 4000h = 0.5 

stl A, *SP (8) ; (AL) --> gain_delta 


ssbx FRCT 
;compute signal power 
ld *SP(9),A 


stlm A, AR2 ;AR2 points to indata2 buffer 


rptz B, # (N-1) ; for (n=0;n<N-1;n++) 
squra *AR2+,B + (B)=(B) +x (n) *x (n) 
sfta B,-6 7 (B)=(1/64) * (B) 
sth B, *SP (0) ; (B) --> gain_pow (32bit) 
stl: B,*SP:(1) ; gain_pow = (N/64) *sigpower 
ld *SP(2),16,A ; (A) = gain_lim = 1/ (64*N) (32bit) 
or *SP(3),A 
sub A,B ;(B) = gain_pow - gain_lim 
be gain2,BLEQ ;if(gain_pow > gain_lim) 
i{ 
; compute gain_amp = sqrt (gain_pow) 


st #4000h,*SP(4) ; initialize gain_amp estimate 
st #4000h,*SP(8) ; initialize gain_delta to 0.5 
mvmm SP,AR2 

mar *+AR2 (8) A AR2 points to gain_delta 

stm #(sqrt_iterations-1),BRC 

rptbd sqrtloopend-1 


squr *SP(4),A ; (A) = yold*yold 
dsub *SP(0),A 7 (A) = yold*yold - x 
sqrtloop 
bed sqrt1,AGT 
ld *SP(4),16,B H (B) = yold 
sub *AR2,15,B ; (B) = yold — (1/2)*delta 
; if(y needs to be larger) 
add *AR2,16,B : (B) = yold + (1/2)*delta 
sqrtl sth B, *SP (4) ; (BH) --> ynew 
ld *AR2,15,A ; (AH) = (1/2)delta 
sth A, *AR2 . (AH) --> delta new 
squr *SP(4),A 
dsub *SP(0),A : (A) = yold*yold - x 
sqrtloopend 
; compute gain_amp = sqrt (64/N) *gain_amp 
ld *SP(7),T A (T) = gain_const 
mpy *SP(4),A ; (A) =sqrt (64/102) *gain_amp 
sth A, *SP (4) . (AH) --> gain_amp 
; compute gain_lev/gain_amp 
ld *SP(5),16,A : (AH) = gain_lev 
rpt #(16-1) 
subc *SP(4),A 7 compute (gain_lev/gain_amp) 
and #0ffffh,A . retain quotient (AL), mask remainder 
sfta A,15 ; shift quotient into high accu 
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sth A, *SP (6) A (AH) --> gain_scale 
; scale indata with (gain_lev/gain_amp) 
ld *SP(9),A , 


stlm A,AR2 , AR2 points to indata 
mvmm SP,AR3 ; 
mar *+AR3 (6) ; AR3 points to gain_scale 
stm #(N-1),BRC ; 
rptb gainl-1 ; for (n=0;n<N;n++) 
mpy *AR2,*AR3,A H (A) = x(n) * (gain_lev/gain_amp) 
sth A, *AR2+ ; (AH) --> x(n) 
gainl Al 


gain2 frame #10 


nop 
popm AR7 
popm AR6 
popm AR1 
popm ST1 
popm STO 
ret 


KKRKKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KK KKKKKKKKKK KKK KKK KKK KKK KKK KKK 


* 


SUBROUTINE: goertzel (w/o OV check approx. 8861 cycles) 


Input: none 
Output: none 
Uses: 


Description: 

This new version of the subroutine operates on a block of data 
102 samples long and computes the following difference equation 
It does thereby not check for overflow and assumes a prescaling 
gain control function to protect from possible overflow. 


vk(n) = 2*coef*vk(n-1) - vk(n-2) + x(n) 


Coefficients are in order starting with location COEF1st 
COEF1st .word XXXX ,coefl 
.word yyyy ;coef2 


Data vk(n-1), vk(n-2) is ordered as follows 
vil5(n-1) 
v15(n-2) 
vl14(n-1) 
v14(n-2) 
v0 (n-1) 
v0 (n-2) <--AR2 


FRCT=1, SXM=1 
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Kkk kkk RRR ARKRKK KKK KKK KKK KKK KKK KKK kkk kkk kkk kkk kkk kkk kkk kkk kkk Kk Kk KK Kk Kk KK Kk KKK 


_goertzel: 
nop 
pshm STO 
pshm STL 
pshm ARI 
pshm AR6 
pshm AR7 
rsbx OVA 
ssbx FRCT 
stlm A,AR5S 
ld *AR5,A 
stlm A,AR1 
ld *AR5(1),A 
add #31,A 
stlm A,AR6 


stm #(N-1),BRC 
goerl rptb goer2-1 
AR6,AR2 
stm #COEF1st,AR3 
ld *AR1+,16,A 


mvmm 


sub *AR2-,16,A,B 
mac *AR2,*AR3,B 
mac *AR2,*AR3+,B 
delay *AR2 
sth B, *AR2- 


sub *AR2-,16,A,B 
mac *AR2,*AR3,B 
mac *AR2,*AR3+,B 
delay *AR2 
sth B, *AR2- 


sub *AR2-,16,A,B 
mac *AR2,*AR3,B 
mac *AR2,*AR3+,B 
delay *AR2 
sth B, *AR2- 


;—---- arguments 


7 (A) = ptr to DTMFCHANNEL 


7;clear overflow 
;set fractional 
;AR5 is used as 


bit for A 
mode 
pointer to struct elements 


; *AR5(0) = 
; *AR5S(1) = 
; *AR5(2) = 
; *AR5(3) = 
; *AR5(4) = 
; *AR5S(5) = 
; *AR5(6) = 
, 

;AR1 points 
; 

; 


;AR6 points 


; for (n=O; n< 
; AR2 poi 
; AR3 poi 
7 (A) = x 
i (B) = - 
i (B) = Cc 
(B) = 2 
; vk (n-1) 
i (B) 

i (B) = = 
i (B) = c 
i (B) = 2 
7 vk (n-1) 
i (B) 

i (B) == 
i (B) = c 
i (B) = 2 
: vk (n-1) 
i (B) 


ptr to indata 

ptr to taps 

ptr to energy template 
digitptr 

digitlast 

detectstat 

err_flags 


to indata buffer 


to end of taps block 


N;nt++) { 
nts to end of taps block 
nts to Coefficients 
(n) 
filter k=0 
vk (n-2) + x(n) 
oef*vk(n-1) - vk(n-2) +x(n) 
*coef*vk(n-1) - vk(n-2) +x(n) 
---> vk(n-2) 
---> vk(n-1) 
filter k=1 
vk (n-2) + x(n) 
oef*vk(n-1) - vk(n-2) +x(n) 
*coef*vk(n-1) - vk(n-2) +x(n) 
---> vk(n-2) 
---> vk(n-1) 
filter k=2 
vk (n-2) + x(n) 
oef*vk(n-1) - vk(n-2) +x(n) 
*coef*vk(n-1) - vk(n-2) +x(n) 
---> vk(n-2) 
-—--> vk(n-1) 
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sub *AR2-,16,A,B ; 
mac *AR2,*AR3,B ; 
mac *AR2,*AR3+,B ; 
delay *AR2 . 
sth B, *AR2- 4 

, 
sub *AR2-,16,A,B 3 
mac *AR2,*AR3,B ; 
mac *AR2,*AR3+,B ; 
delay *AR2 ; 
sth B, *AR2-— 7 

r 
sub *AR2-,16,A,B ; 
mac *AR2,*AR3,B ; 
mac *AR2,*AR3+,B 7 
delay *AR2 ; 
sth B, *AR2- ’ 

, 
sub *AR2-,16,A,B i. 
mac *AR2,*AR3,B ; 
mac *AR2,*AR3+,B ; 
delay *AR2 ; 
sth B, *AR2- ; 

, 
sub *AR2-,16,A,B ; 
mac *AR2,*AR3,B ; 
mac *AR2,*AR3+,B ; 
delay *AR2 7 
sth B, *AR2- ; 

, 
sub *AR2-,16,A,B i 
mac *AR2,*AR3,B ; 
mac *AR2,*AR3+,B ; 
delay *AR2 iF 
sth B, *AR2- ; 

, 
sub *AR2-,16,A,B - 
mac *AR2,*AR3,B 7 
mac *AR2, *AR3+,B i: 
delay *AR2 ; 
sth B, *AR2- ; 

, 
sub *AR2-,16,A,B ; 
mac *AR2,*AR3,B f 
mac *AR2,*AR3+,B ; 
delay *AR2 ; 
sth B, *AR2- 7 

, 
sub *AR2-,16,A,B : 
mac *AR2,*AR3,B ’ 
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Sen Ss= filter k=3 ----------------- 
(B) = -vk(n-2) + x(n) 

(B) = coef*vk(n-1) - vk(n-2) +x(n) 
(B) = 2*coef*vk(n-1) - vk(n-2) +x(n) 
vk(n-1) ---> vk(n-2) 

(B) -—--> vk(n-1) 

=== filter k=4 ----------------- 
(B) = -vk(n-2) + x(n) 

(B) = coef*vk(n-1) - vk(n-2) +x(n) 
(B) = 2*coef*vk(n-1) - vk(n-2) +x(n) 
vk(n-1) ---> vk(n-2) 

(B) —--> vk(n-1) 

=SSs=> filter k=5 ----------------- 
(B) = -vk(n-2) + x(n) 

(B) = coef*vk(n-1) - vk(n-2) +x(n) 
(B) = 2*coef*vk(n-1) - vk(n-2) +x(n) 
vk(n-1) ---> vk(n-2) 

(B) —--> vk(n-1) 

SSSSS> filter k=6 ----------------- 
(B) = -vk(n-2) + x(n) 

(B) = coef*vk(n-1) - vk(n-2) +x(n) 
(B) = 2*coef*vk(n-1) - vk(n-2) +x(n) 
vk(n-1) ---> vk(n-2) 

(B) —--> vk(n-1) 

<= filter k=7 ----------------- 
(B) = -vk(n-2) + x(n) 

(B) = coef*vk(n-1) - vk(n-2) +x(n) 
(B) = 2*coef*vk(n-1) - vk(n-2) +x(n) 
vk(n-1) ---> vk (n-2) 

(B) —--> vk(n-1) 

SSSSS= filter k=8 ----------------- 
(B) = -vk(n-2) + x(n) 

(B) = coef*vk(n-1) - vk(n-2) +x(n) 
(B) = 2*coef*vk(n-1) - vk(n-2) +x(n) 
vk (n-1) ---> vk (n-2) 

(B) —--> vk(n-1) 

tt filter k=9 ----------------- 
(B) = -vk(n-2) + x(n) 

(B) = coef*vk(n-1) - vk(n-2) +x(n) 
(B) = 2*coef*vk(n-1) - vk(n-2) +x(n) 
vk(n-1) ---> vk (n-2) 

(B) —--> vk(n-1) 

=====> filter k=10 ----------------- 
(B) = -vk(n-2) + x(n) 

(B) = coef*vk(n-1) - vk(n-2) +x(n) 
(B) = 2*coef*vk(n-1) - vk(n-2) +x(n) 
vk (n-1) ---> vk (n-2) 

(B) ---> vk(n-1) 

SSSSS> filter k=11 ----------------- 
(B) = -vk(n-2) + x(n) 

(B) = coef*vk(n-1) - vk(n-2) +x(n) 


mac *AR2,*AR3+,B ; (B) = 
delay *AR2 ; vk (n-1) 
sth B, *AR2- H (B) 
e.  “Seteee 
sub *AR2-,16,A,B ; (B) = 
mac *AR2,*AR3,B A (B) = 
mac *AR2,*AR3+,B ; (B) = 
delay *AR2 7 vk (n-1) 
sth B, *AR2- ; (B) 
eee 
sub *AR2-,16,A,B H (B) = 
mac *AR2,*AR3,B ; (B) = 
mac *AR2, *AR3+,B i (B) = 
delay *AR2 ; vk (n-1) 
sth B, *AR2- ; (B) 
See 
sub *AR2-,16,A,B ; (B) = 
mac *AR2,*AR3,B : (B) = 
mac *AR2,*AR3+,B H (B) = 
delay *AR2 ; vk (n-1) 
sth B, *AR2- ; (B) 
Be ee 
sub *AR2-,16,A,B ; (B) = 
mac *AR2,*AR3,B ; (B) = 
mac *AR2,*AR3+,B ; (B) = 
delay *AR2 ; vk (n-1) 
sth B, *AR2- ? (B) 
goer2 i} 
goer3 popm AR7 

popm AR6 

popm AR1 

popm STl 

popm STO 

ret 


2*coef*vk (n-1) 
---> vk(n-2) 
-—--> vk(n-1) 

filter k=12 

-vk(n-2) + x(n) 

coef*vk (n-1) 

2*coef*vk (n-1) 
-—--> vk(n-2) 
-—--> vk(n-1) 

filter k=13 

-vk(n-2) + x(n) 

coef*vk (n-1) 

2*coef*vk (n-1) 
---> vk(n-2) 
---> vk(n-1) 

filter k=14 

-vk(n-2) + x(n) 

coef*vk (n-1) 

2*coef*vk (n-1) 
---> vk(n-2) 
---> vk(n-1) 

filter k=15 

-vk(n-2) + x(n) 

coef*vk (n-1) 

2*coef*vk (n-1) 

---> vk(n-2) 

-—--> vk(n-1) 


— vk (n-2) 


— vk (n-2) 
— vk (n-2) 


— vk (n-2) 
— vk (n-2) 


— vk (n-2) 
— vk (n-2) 


— vk (n-2) 
— vk (n-2) 
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+x (n) 


+x (n) 
+x (n) 


+x (n) 
+x (n) 


+x (n) 
+x (n) 


+x (n) 
+x (n) 


Kkk kkk kK KR ARKARK KKK KKK KKK KKK Kk kK KKK KKK kkk kkk kk kkk kkk kkk kK kkk kK KKK Kk KKK KK kk KK KKK 


* 


* 


* 


* 


* 


SUBROUTINE: dtmf_checks 


Input: none 

Output: A (status if valid digit) 
Uses: 

Description: 


— energy computation 

-— initialize filter taps to zero 
— maximum search 

- DIMF status check 

- signal strength check 

- twist check 

- relative peak check 

- 2nd harmonic check 

- digit validation 


(approx 788 cycles) 
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* 


KKKKKKKKKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KK KKKKKKKKKKK KKK KKK 


_dtmfchecks: 
nop 
pshm STO 
pshm ST1 
pshm AR1 
pshm AR6 
pshm AR7 


nop 
frame #0 ;---- local variables ----- 
; none 
{>> argunents: --—---=-=+= 
7 (A) = ptr to dtmfchannel struct 
ssbx FRCT 
stlm A,ARS ;AR5 is used as ptr to dtmfchannel struct 
7 *AR5(0) = ptr to indata 
, *AR5(1) = ptr to taps 
; *AR5(2) = ptr to energy template 
; *AR5(3) = digitptr 
; *AR5(4) = digitlast 
7 *AR5(5) = detectstat 
; *AR5(6) = err_flags 
j -------------------------------------~---------~-~----~-~--~-~------~---- 
; Energy computation 
; 
7 Y(N)y*(N) = vk(N) *vk(N) - 2*coef*vk (N) *vk(N-1) + vk(N-1) *vk (N-1) 


; Coefficients are in pmem in order starting with location COEFI1st 
; COEF1st .word RHRH ;coef0 <--AR3 
; -word  yyyy ,coefl 


; Data vk(n-1), vk(n-2) is ordered as follows in dmem 
; vi5(n-1) 


H v1l5(n-2) 
H v14(n-1) 
s vl14(n-2) 
: v0 (n-1) 


iH v0 (n-2) <--AR2 


; FRCT = 1, OVLY =1 
; AR2: points to v1 (N-1) 
; AR3: points to COEF1st 


ener ld *AR5(1),A 


add #31,A 
stlm A,AR2 ;AR2 points to end of filter taps block 
stm #COEF1st,AR3 ;AR3 points to beg of coefficient table 
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ld *AR5(2),A 


stlm A,AR4 ;AR4 points to energy template 

stm #(16-1),BRC 
rptb enerl-1l ; for (k=0;k<16; k++) { 

ld *AR2,16,A ; (A) = vk(N-1), A(32:16) loaded 

mpya *AR2— 7 (B) = vk(N-1)*vk(N-1), (T) = vk(N-1) 

mpy *AR2,A 7 (A) = vk(N) *vk (N-1) 

ld *AR3+,T > (T) = coef 

mpya A ; (A) = coef*vk (N) *vk (N-1) 

sub A,1,B ; (B) = -2*coef*vk (N) *vk(N-1) + vk (N-1) *vk (N-1) 
ld *AR2,T ; (T) = vk(N) 

mac *AR2-,B 7 (B) = vk(N)*vk(N) - 2*coef*vk (N) *vk(N-1) + vk (N-1) *vk (N-1) 
sth B, *AR4+ ; store result E[k] into energy data block 

enerl i} 


ld *AR5(1),A 


stlm A,AR2 ;AR2 points to beginning of filter taps 
ld #0,A 
rpt #(32-1) ; for (i=0;1i<32; i++) 
stl A, *AR2+ ; initialize filter taps for next round 
BREAK nop ;<--- breakpoint for testing 


; Alternative Maximum search 
7 result: AR3 points to max row energy ROWMAX 
7 AR4 points to max col energy COLMAX 
max_search: 

ld *AR5(2),A 


stlm A,AR2 ;AR2 points to ENERGY table (row energies) 
mvmm AR2,AR3 ;AR3 points to ENERGY table 
stm #(4-1),BRC 
rptb maxl-1 ; for (k=0;k<4; k++) 
ld *AR3,A ; (A) =oldmax 
sub *AR2,A : (A) =oldmax-E[k] 
bce maxl11,AGT iH if (oldmax < E[k]) 
mvmm AR2,AR3 ; AR3 points to newmax 
maxll mar *AR2+ A increment pointer 
maxl ;result: AR3 points to ROWMAX 
; AR2 points to column energies 
mvmm AR2,AR4 ;AR4 points to ENERGY table (column energies) 
stm #(4-1),BRC 
rptb max2-1 ; for (k=0;k<4; k++) 
ld *AR4,A : (A) =oldmax 
sub *AR2,A : (A) =oldmax-E[k] 
be max21,AGT : if(oldmax < E[k]) 
mymm AR2,AR4 ; AR4 points to newmax 
max21 mar *AR2+ . increment pointer 
max2 ;result: AR4 points to COLMAX 
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; AR2 points to column energies 


; DIMF Status Check: 


; if (detectstat == 0) 

; if ( max (ROWMAX,COLMAX) < THR_PAU ) 
; detectstat = 1, detector enabled 
; terminate all checks 

7 else 

; continue 


; AR3 points to ROWMAX 

; AR4 points to COLMAX 

stat_check 
andm #0080h, *AR5(6) ;clear error flags, except OVA 
bitf *AR5(5),#0001h ; 


be sig_check,tc ;if(detectstat == 1) continue with sig_check 
ld *AR3,A ; 
ld *AR4,B ;else 
max A H (A) =max (ROWMAX, COLMAX) 
sub #THR_PAU,A ; (A) =max (ROWMAX, COLMAX) - THR_PAU 
be statl,ageq ;  if((A)<0) 
st #0001h, *AR5 (5) ; detectstat = 1 
statl b ende H terminate all checks 


; Signal strength check: 
; if (ROWMAX+COLMAX > THR_SIG) continue 
; else terminate all checks 


; AR3 points to ROWMAX 
; AR4 points to COLMAX 


sig_check: 


ld *AR3,A 

add *AR4,A ; (A) = ROWMAX+COLMAX 

sub #THR_SIG1,A,B ; 

bce ende, BLEQ ; if (ROWMAX+COLMAX <= THR_SIG1) 


: signal is not possible tone 
; terminate checks 
sub #THR_SIG2,A,B ; 
be err_sig, BLEQ ;if (ROWMAX+COLMAX <= THR_SIG2) 
i; signal is possible tone 
4 however does not meet THR_SIG2 
: branch to err_sig 
b twist_check 


err_sig: 7**** record signal strength error **** 
orm #0001h,*AR5(6) ; set signal strength error flag (b0) 
b ende ; terminate checks 


; Twist check: 
; if (ROWMAX > COLMAX) 
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; check reverse twist: 

7 if ((COLMAX/ROWMAX) > THR_REVTWI) continue 
} else terminate all checks 

; else 

; check standard twist: 

7 if ((ROWMAX/COLMAX) > THR_STDTWI) continue 
; else terminate all check 


; AR3 points to ROWMAX 
; AR4 points to COLMAX 


twist_check: 
ld *AR3,A 
sub *AR4,A 
be twist2,ALEQ 


;if (ROWMAX > COLMAX) 


twistl ld *AR4,16,A Fd (AH) =COLMAX 
rpt #(16-1) ; 
subc *AR3,A j compute (COLMAX/ROWMAX) 
and #0ffffh,A ; retain quotient (AL), mask remainder 
sfta A,15 ; shift quotient into AH 


sub #THR_REVTWI,16,A; 


be err_revtwi,ALEQ 7 if ((COLMAX/ROWMAX) <= THR_REVTWI) 


H branch to err_revtwi 
b rel_check 
7else 
twist2 ld *AR3,16,A : (AH) =ROWMAX 
rpt #(16-1) i 
subc *AR4,A j compute (ROWMAX/COLMAX) 
and #0ffffh,A ; retain quotient (AL), mask remainder 
sfta A,15 : shift quotient into AH 


sub #THR_STDTWI,16,A_ ; 


be err_stdtwi,ALEQ ; if ((COLMAX/ROWMAX) <= THR_STDTWI) 


; branch to err_stdtwi 

b  rel_check 

err_revtwi: ;**** record reverse twist error **** 
orm #0002h, *ARS (6) ; set reverse twist error flag (bl) 
bende ; terminate checks 

err_stdtwi: ;**** record standard twist error **** 
orm #0004h, *AR5 (6) ; set standard twist error flag (b2) 
bende ; terminate checks 

j ----------------------------------------------------~-~---~-~-~--~--~--- 

; Relative peak check: 

; 

; find rel row peak -—-> RELPEAK 

; if (RELPEAK/ROWMAX < THR_ROWREL) continue 

; else terminate all checks 

, 

; find rel col peak --> RELPEAK 

; if (RELPEAK/COLMAX < THR_ROWREL) continue 

; else terminate all checks 
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ts 


AR3 points to ROWMAX 
AR4 points to COLMAX 


rel_check: poccoee relative row peak ratio check ------- 
ld *AR5(2),A ; 
stlm A, ARO ;ARO points to ENERGY table (row frequencies) 
1d #0,B ; (BH) =0 
stm #(4-1),BRC ; 
rptb rel12-1 ; for (k=0;k<4; k++) 
ld *ARO,16,A ; (AH) =E [k] 
cmpr 0,AR3 ; 
be rell11,TC ; if (ARO==AR3) 
max B 7 (B) =newmax=max (oldmax, E[k] ) 
relll mar *ARO+ ; 
rell2 rpt #(16-1) : 
subc *AR3,B ;compute (RELPEAK/ROWMAX) 
and #0ffffh,B ;retain quotient (BL), mask remainder 
sfta B,15 ;shift quotient into BH 


sub #THR_ROWREL,16,B ; 
be err_rowrel, BGEQ ; if ((RELPEAK/ROWMAX) >= THR_ROWREL) 


; branch to err_rowrel 
j=-=-=- relative column peak ratio check -—--- 
rel2 ld #0,B ; (BH) =0 
stm #(4-1),BRC ; 
rptb rel22-1 ; for (k=4;k<8; k++) 
ld *ARO,16,A ; (AH) =E([k] 
cmpr 0,AR4 ; 
be rel21,TC ; if (ARO==AR4) 
max B 7 (B) =newmax=max (oldmax, E[k] ) 
rel21 mar *ARO+ ; 
rel22 rpt #(16-1) : 
subc *AR4,B ;compute (RELPEAK/COLMAX) 
and #0ffffh,B ;retain quotient (BL), mask remainder 
sfta B,15 ;shift quotient into BH 


sub #THR_COLREL,16,B ; 
be err_colrel,BGEQ ; if ((RELPEAK/COLMAX) >= THR_COLREL) 


; branch to err_colrel 
b sec_check 
err_rowrel: 7**** record row’s rel peak error **** 
orm #0008h, *AR5 (6) ; set row’s rel peak error flag (b3) 
b ende ; terminate checks 
err_colrel 7**** record col’s rel peak error **** 
orm #0010h, *AR5 (6) ; set col’s rel peak error flag (b4) 
b ende ; terminate checks 
j ------------------------------------~-----------~----~-~----~-----~--~--- 
; Second Harmonics check 
> if( (ROW2nd/ROWMAX) < THR_ROW2nd ) continue 


else terminate all checks 
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7 if( (COL2nd/COLMAX) < THR_COL2nd ) 
} else terminate all checks 


7 AR3 points to ROWMAX 
; AR4 points to COLMAX 
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continue 


sec_check: 
stm #8h, ARO 


mar *AR3+0 


ld 


*AR3-0,16,A 


sub *AR3,16,A,B 


be 


err_row2nd, BGEQ 


rpt #(16-1) 
subc *AR3,A 

and #0ffffh,A 
sfta A,15 


;offset into 2nd harmonics energies 

== =SS2 2nd harmonics check for rows -—------— 
;AR3 points to corresponding ROW2nd 

; (AH) =ROW2nd, AR3 points to ROWMAX 

; (BH) =ROW2nd-ROWMAX 

7; if (ROW2nd<ROWMAX) 

7 

; compute (ROW2nd/ROWMAX) 

7 retain quotient (AL), mask remainder 

; shift quotient into AH 


sub #THR_ROW2nd,16,A_ j; 
be err_row2nd,AGEQ ; 


if ( (ROW2nd/ROWMAX) 


>= THR_ROW2nd) 


branch to err_row2nd 


sec2 
ld 


mar *AR4+0 
*AR4-0,16,A 


sub *AR4,16,A,B 


be 


err_col2nd, BGEQ 


rpt #(16-1) 
subc *AR4,A 

and #0ffffh,A 
sfta A,15 


sub #THR_COL2nd,16,A 


be 


b 


err_col2nd, AGEQ 


digit_map 


err_row2nd: 
orm #0020h, *AR5 (6) 


b 


ende 


err_col2nd: 
orm #0040h, *AR5 (6) 


b 


ende 


; 2nd harmonics check for columns ----— 
;AR4 points to corresponding COL2nd 

; (AH) =COL2nd, AR4 points to COLMAX 

; (BH) =COL2nd—-COLMAX 

; if (COL2nd<COLMAX) 

r 

; compute (COL2nd/COLMAX) 


; retain quotient (AL), mask remainder 


H shift quotient into AH 


’ if ((COL2nd/COLMAX) >= THR_COL2nd) 
7 branch to err_col2nd 


7;**** record row’s 2nd harm error **** 
; set row’s 2nd harm error flag (b5) 

; terminate checks 

7 **** record col’s 2nd harm error **** 
; set col’s 2nd harm error flag (b6) 

; terminate checks 


; Map detected tone pair to digit 


; AR3 points to ROWMAX 
; AR4 points to COLMAX 


find row and col numbers -—--~-~--~-~-~--~--~--~--~-----------— 


digit_map: 
j---------- 
ld *AR5(2),A 
stlm A, ARO ;load offset 
nop 
nop 
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mar *AR3-0 ;AR3 contains row# 
mar *+ARO (4) 
mar *AR4-0 ;AR4 contains col# 
=== — map row and column numbers to actual digit ------------- 
ld *(AR3),8,A ;load row nbr in upper byte of AL 
add * (AR4),A ;load column nbr in lower byte of AL 
stm #KEYS, AR2 ;AR2 points to key mapping table 
stm #0,AR3 ;AR3 is used as counter 
stm #(16-1),BRC 
rptb digit3-1 ; for (k=0;k<16; k++) { 
sub *AR2+,A,B 
be digit31,bneq ; if (A == KEYS[k]) insert delay 
b digit4 ; branch to ende 
digit31 mar *AR3+ : increment counter 
digit3 i} 


digit4 ;now AR3 contains decoded digit 
ld *AR5(3),A 
stim A, AR4 7AR4 used as digitptr 
ldm AR3,A ;load current digit 
sub *AR5(4),A,B ;compare with last digit 
be digit5,bneq ;if current digit and last digit are the same 
stl A, *AR4+ ; store VALID DIGIT into DIGITS 
1ldm AR4,A ; 
stl A, *AR5 (3) ; update pointer 
st #0, *AR5 (5) ; disable detector and wait for pause 
digit5 ldm AR3,A ; store decoded digit 
stl A, *AR5 (4) ; into digitlast 
ende frame #0 


popm AR7 
popm AR6 
popm AR1 
popm ST1 
popm STO 
nop 

ret 
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COPYRIGHT TEXAS INSTRUMENTS, 


INC. 


1996 


* 
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m (C) 

im Program Name: 

- File Name: 

ms File Description: 


* Author: 


* Date: 10/15/96 
* Revision: 2.0 
* Latest working date: 


Gunter Schmer 


DIMF tone decoder 
tables.inc 


tables for dtmfsub.asm 


10/15/96 


KKK KK KKK KK KKK KKK KKK KKK KKK KKK KKK KKK KEKE KEKE KEKE KEKE KKK KKKKK KKK KKK KKK KKK 


KKK KKK KKK KKK KEK KKK KEKE KKKKKKK KKK KKK KK KK KKK 


KkK* 


Constants 


Kak kK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK Kk KK 


sqrt_iterations 


N «Seu 

THR_SIG1 -set 
THR_SIG2 -set 
THR_PAU -set 
THR_STDTWI .set 
THR_REVTWI .set 
THR_ROWREL .set 
THR_ROW2nd .set 
THR_COLREL .set 
THR_COL2nd .set 


102 


1900*32768/10000 
2150*32768/10000 
458*32768/10000 
3162*32768/10000 
1259*32768/10000 
2000*32768/10000 
2000*32768/10000 
2000*32768/10000 
2000*32768/10000 
-set 16 


7;max count 


KKK KKK KKK KKK KKK KK KEK KKK KKKK KKK KKK KK KKK KK 


mem Tables 
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;threshol 
;threshol 
;threshol 
;threshol 
;threshol 
;threshol 
;threshol 
;threshol 
; threshold 
;iterations for sqrt 


ld 
ld 
ld 
ld 
ld 
ld 
ld 
ld 


£or 
for 
for 
for 
for 
for 
for 
for 
for 
for 


* 


* 


DFT loops 
possible tone 


definite tone 
pause energy 
standard twist (-5dB) 
reverse twist (-9dB) 


row’s relative peak ratio 
row’s 2nd harmonic ratio 
col’s relative peak ratio 
col’s 2nd harmonic ratio 


in gain control function 


Kkhkh kkk kK RAR ARK K KKK KKK KKK KKK Kk kkk kk kkk kkk Kk Kk KK KK KK KKK KKK KKK 


* The following table holds the coefficients for 


* the 16 Goertzel filters, 


* 


- -word 

* 

* with WN = 102 
“ k= 9, 
* k = 18, 


* 


cos (2pi*k/N) *32768 


(lst and 2nd Harmonics) 


10, 
20, 


11, 
22, 


13, 
24, 


15, 
31, 


17, 
34, 


19, 
38, 


21 
42 


computed as follows 


(dhs 
(25 


; coefficient 


Harm) * 


Harm) * 
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"tbhl_coeft” 


-sect 
COEF1st .word 27980 
-word 26956 
-word 25701 
-word 24219 
-word 19073 
-word 16325 
-word 13085 
-word 9315 
COEF2nd .word 14739 
-word 11414 


, 


’ 


1st Harmonics 


2nd Harmonics 
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-word 7549 


-word 3032 

sword -10565 
-word -16503 
-word -22318 
-word -27472 
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* Dis 
* 


* 


git mapping table 
Low byte row number 


High byte: column number 


* 


* 


* 
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KEYS 


-sect "tbl_keys” 
-word 0301h 
-word 0000h 
-word 0001h 
-word 0002h 
-word 0100h 
-word 0101h 
-word 0102h 
-word 0200h 
-word 0201h 
-word 0202h 
-word 0003h 
-word 0103h 
-word 0203h 
-word 0303h 
-word 0300h 
-word 0302h 


, 


’ 


ro: 
my 
ror 
13° 
rar 
or 
"6! 
yr 
rer 
ror 
ar 
BI 
rc 
pr 
Er 
Br 


a 


rye 
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* Test Buffer for encoding results 
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DATA 


D-18 


-usect "buffer”,8000h ;length 32k words 
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* 


